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TECHNICAL MEMORANDUM 


A “KANE’S DYNAMICS” MODEL 
FOR THE ACTIVE RACK ISOLATION SYSTEM 

PART TWO: NONLINEAR MODEL DEVELOPMENT, 
VERIFICATION, AND SIMPLIFICATION 


1. INTRODUCTION 


In part one of this series, NASA/TM— 2001-211063, a high-fidelity, linearized, analytical model 
of the active rack isolation system (ARIS) was derived using traditional, hand-calculation methods. 1 
The model was developed directly, using Kane’s method, without intermediate development of a full 
nonlinear model. 2 ’ 3 This Technical Memorandum (TM), part two, presents four computer-based, num- 
erical models of ARIS, one of which is purely kinematical, and the remaining three, dynamical. These 
numerical models are used collectively to verify the linearized analytical model developed in part one. 

The kinematical model uses high-fidelity computer-aided design (CAD) models and contains 
17 rigid bodies and 48 degrees of freedom. It uses specified positions of actuator number 1, and 
ARIS geometry, to predict the configurations of the remaining seven actuators. 

The three dynamical models are a simple truth model, a full nonlinear rigid-body model without 
umbilicals, and a corresponding linearized rigid-body model. The first dynamical model represents the 
system dynamics of an international standard payload rack (ISPR) in zero gravity with characteristic 
mass properties but without any actuators; it is used as a truth model, providing a benchmark for com- 
parison purposes. The second dynamical model represents the full nonlinear system dynamics of an 
ARIS-outfitted ISPR. The actuator spring constants in the nonlinear model are first reduced, from nomi- 
nal to negligible levels, for comparisons against the truth model using a set of test inputs. The actuator 
spring constants are then increased to moderate levels, to provide a more realistic representation of the 
system dynamics. The third dynamical model is a numerical model that has been linearized from the 
nonlinear model; it is a numerical representation of the linearized analytical model presented in part one 
The input responses of this linearized model are compared to those of the nonlinear model. These com- 
parisons confirm that the former, like the latter, accurately represents the physics of the system for ISPR 
motions within operational limits. 

Once the three dynamical models have been run and compared, the kinematical model is 
used for visualization. The kinematical model drives a CAD representation of the ARIS config 
uration for specified times during the simulation. The result is a clear, visual representation of 
dynamical-model predictions of the configuration for the entire system at the various, arbitrary 
stages during the simulation. 


2. DESCRIPTION OF SOFTWARE 


The development and verification of the dynamical models were accomplished using two com- 
mercial-off-the-shelf (COTS) software tools: OnLine Dynamics’ Autolev™ and DELMIA’s® (formerly 
DENEB Robotics’) ENVISION . 3 - 5 

ENVISION is a robotics software package developed for the automotive industry and employs 
three-dimensional CAD models to facilitate both forward and inverse kinematics analyses. Once CAD 
models of all of the rigid bodies of ARIS have been translated into ENVISION, coordinate systems 
representing the attachment locations are added. Then, in a parent-child fashion, from the CAD models, 
each ARIS actuator is built, and an algorithm for an inverse kinematics solution is applied to each. Next, 
all of the actuators, the International Space Station ( ISS ), and ISPR are appropriately combined, and the 
system is activated in a common CAD environment. 

At this point, a user can specify the configuration of a single input actuator, and the remain 
ing seven actuators are forced to their unique configurations as determined by ENVISION’s inverse 
kinematics algorithm and the constraints of the system. The configuration information for each of 
the actuators can then be compared to the corresponding configuration indicated by an Autolev 
simulation. 

Autolev is a DOS-based interpreter that is designed, in general, to solve vector-based mathemat- 
ical problems and, in particular, solve dynamics problems using Kane’s method. Because Autolev is an 
interpreter, all of the commands entered for a particular session are stored in memory. A user can easily 
develop a text file containing the commands for a particular session. These files can be read by Autolev 
as if the commands had been entered manually. 

The output from Autolev is a C program. This program when compiled, linked, and run, creates 
a DOS-executable file that, in turn, creates a data file of the specified output for a given simulation run. 
These automatically generated files are not included in this TM since they are very long and not easy 
to follow. Most of the descriptions in this TM illustrate the application of Autolev to the ARIS problem. 
There is a close correlation between the methods for the development of the closed-form analytical 
ARIS model and the Autolev model because both use Kane’s method. 
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3. TRUTH MODEL 


The first Autolev model represents the system dynamics of an ISPR in zero gravity with charac- 
teristic mass properties but without any actuators. Developed in Autolev, this model tracks the response 
in translation and rotation of the ISPR center of mass, as test forces are applied at various points on 
ISPR. The coordinates used in this model map the translation along and rotation about 
the x-, y-, and z-axes. 

This model serves as a baseline for the comparison of other, more complex models, primarily 
because hardware validation would require resources beyond those currently available. One positive 
aspect of this model is that it is very simple so that its results are easily reproduced using various tech- 
niques. A limitation of this model is that it is useful only when the dynamical effects of the actuators 
are negligible. This requirement means, first, that the torsional springs at the lower and upper stingers 
and cross flexures must have a negligible stiffness, and second, the actuators must not reach full 
extension. 

The actuator spring constants are easily changed to negligible levels. This is accomplished 
by specifying very small values for the spring constants in the code. The full-extension constraints 
only come into play when the actuators would otherwise extend beyond their physical limitations. 

In the actual hardware, the ISPR displacements from the home position are limited by a snubber sys- 
tem. However, in the Autolev model, there is no such system, so the motion must be limited by the dis- 
turbance input. 
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4. NONLINEAR ACTIVE RACK ISOLATION SYSTEM MODEL 


Kane’s method was used to develop both the linearized analytical model, presented in part one 
of this series, and the corresponding Autolev model discussed in this TM. 1-3 Most of the developmen- 
tal steps outlined in part one are implemented specifically by corresponding commands in Autolev, 
and in a few cases, Autolev automates several steps with a single command. A copy of the text file 
that represents the nonlinear ARIS model is given in A High Fidelity Model for the Active Rack Isola- 
tion System f 


4.1 Rotation Matrices 

Using Autolev notation, let the Ai j > coordinate system rotate, relative to the Fi j > coordi- 
nate system, through positive angle Q ( l + 6n) , (n=0 , 7 ) about the Fi3> axis. Similarly, let the 
orientation of the Ai j > coordinate system, relative to the Pi j > coordinate system, be described by 
consecutive positive rotations : Q ( 2 + 6 n ) about the P i 1 > axis and Q ( 3 + 6 n ) about the moved 3 
axis. Also, let the orientation of the Pi j > coordinate system, relative to the Si j > coordinate system, 
be described by consecutive positive rotations: Q ( 4 + 6 n ) about the S i 3 > axis, Q ( 5 + 6 n ) about the 
moved 2 axis, and Q ( 6+6n) about the moved 1 axis. 

As entered into Autolev, the rotation matrices among the several coordinate systems for the 
first actuator assembly are 


SIMPROT (Sa,Ra, 3, Q4) , 

( 1 ) 

SIMPROT (Sa,Ra, 2, Q5) , 

( 2 ) 

SIMPROT (Qa, Pa, 2, Q6) , 

( 3 ) 

SIMPROT (Pa, Ta, 2, Q2) , 

( 4 ) 

SIMPROT (Ta,Aa, 2, Q3) , 

( 5 ) 


SIMPROT (Fa, Aa, 2 , Q1 ) . (6) 

Define the rotation matrix between each of the eight, stator- fixed coordinate systems, Si j >, and the 
single, stator-fixed, reference coordinate system, S j >. Measurements from the ENVISION model 
for the first actuator show that 

Sa_S=[0, -0. 7071, 0.7071;0, 0.7071, 0.7071;-1, 0, 0] , ( 7 ) 

where the commas separate elements within a row and the semicolons separate the rows. 


4 


Next, define the rotation matrix between each of the eight, flotor-fixed coordinate systems, 

Fi j >, and the single, flotor-fixed, reference coordinate system, Fj >. For the first actuator (from 
the ENVISION model), 

F_Fa=[0, 0,1;0. 86603,-0. 5, 0;0. 5, 0.8663,0] . (8) 

Note that the nomenclature, F_Fa, represents the rotation from the Fa> system, fixed on the 
ISPR at the cross flexure to the F> system, fixed on the flotor at the center of mass. 

4.2 Generalized Coordinates 

The 48 angles, Q j ( j =1 , 4 8 ) , are the generalized coordinates of the system. For the zth 
actuator, the six associated generalized coordinates are as follows: 

• Q ( 1 + 6n ) ( n=0 , ... , 7 ) is the angle at the cross flexure of the ith actuator. 

• Q ( 2 + 6 n ) and Q ( 3 + 6 n ) are the angles at the upper stinger. 

• Q ( 4 + 6n) , Q (5 + 6n) , and Q ( 6+6n) are the angles at the lower stinger. 

4.3 Generalized Speeds 

Define generalized speeds U j ( j =1...4 8 ) for the system as the time rates of 
change of the independent generalized coordinates of S in the inertial reference frame: 

Uj =Qj ' ( j=l . . . 6) . (9) 

4.4 Angular Velocities of the Reference Frames and Rigid Bodies 

Designate the reference frames corresponding to the stator, the zth push rod, the zth arm, and the 
flotor by the symbols Sj>, Pij>,Aij>, and Fj >, respectively. Let Si j > and Fi j > represent, respec- 
tively, the coordinate systems in S j > and Fj > located at the respective points Oi and Fli. 

Two intermediate reference frames were introduced previously to permit describing the angular 
velocity of each pushrod relative to the stator. Designate those intermediate frames corresponding to the 
ith actuator assembly by Ri j > and Qi j >. Another intermediate reference frame was previously intro- 
duced between frames Pi j > and Ai j >; designate this by Ti j >. 

With the expression for angular velocity denoted W_Aa_Fa> corresponding to the angular 
velocity of Aa j > relative to Fa j >, the following expressions represent the angular velocities 
through the first actuator in Autolev: 


1 

00 

1 

2 

V 
II 

o 

V 

(10) 
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1 

00 

0> 

1 

00 

V 
II 

o 
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(ID 

W Ra Sa>=U4*Sa3> 

(12) 
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( 13 ) 


and 


W_Qa_Sa>=U4 *Sa3>+U5*Ra2> , 
W_Pa_Sa>=U4*Sa3>+U5*Ra2>+U6*Qal> , ( 14 ) 

w _ Ta _ Sa>=u4 * Sa3>+u5 * Ra2>+u6 *Q al>+u2 * Pal> ’ (I 5 ) 

W_Aa_Sa>=U4 *Sa3>+U5*Ra2>+U6*Qal>+U2*Pal>+U3*Ta3> , ( 16 ) 

W_Fa_Sa>=U4*Sa3>+U5*Ra2>+U6*Qal> 

+U2*Pal>+U3*Ta3>-Ul*Fa> , ( 17 ) 

W_Aa_Pa>=U2*Pa2>+U3*Ta3> , ( 18 ) 

w _ Fa _ Aa>=_ul * Fa3> ’ ( 19 ) 

W_F_Fa>=0> . ( 20 ) 


4.5 Angular Accelerations of the Reference Frames and Rigid Bodies 

The designation ALF_X_Y> represents the angular acceleration of the coordinate system or 
rigid body, X, with respect to the coordinate system or rigid body, Y. They are calculated as the time 


derivatives of the angular velocities above, using the Autolev commands 

ALF_S_N>=0> , ( 21 ) 

AL F_S a_S >= 0 > , ( 22 ) 

ALF_ R a_Sa>=DT (W_ R a_Sa>, Sa) , ( 23 ) 

ALF_Qa_Sa>=DT (W_Qa_Sa>, Sa) , ( 24 ) 

ALF_Pa_Sa>=DT (W_Pa_Sa>, Sa) , ( 25 ) 

ALF_Ta_Sa>=DT (W_Ta_Sa>, Sa) , ( 26 ) 

ALF_Aa_Sa>=DT (W_Aa_Sa>, Sa) , ( 27 ) 

ALF_Fa_Sa>=DT (W_Fa_Sa>, Sa) , ( 28 ) 

ALF_Aa_Pa>=DT (W_Aa_Pa>, Sa) , ( 29 ) 

ALF_Fa_Aa>=DT (W_Fa_Aa>, Sa) , ( 30 ) 
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and 


ALF F Fa>=0> . 


( 31 ) 


4.6 Velocities of the Centers of Mass for the Rigid Bodies 


Represent by P_B_C> the position vector from arbitrary point B to arbitrary point C. Define 
the following position vectors for the first actuator (and similarly for the remaining seven), as measured 
from the ENVISION model: 


and 


p _ 0a _ A2a>=1 • 375* pa2> , ( 32 ) 

p _ 0a _ Pao>= ° • 6875 * Pa2> > ( 33 ) 

p _ A2a _ Fla>=_0 • 25 * Aal> ’ ( 34 ) 

p _ A2a _ Aao>=0 • 2375* Aal> +0 . 079 1 7* Aa2> , ( 35 ) 

P_Fl a _Fo>=2 . 122*Fal>+l . 578*Fa2>-l . 494*Fa3> . ( 36 ) 


The first derivatives of the appropriate position vectors, under the stated assumptions, yield expres- 
sions for the velocities of the centers of mass for the 17 rigid bodies. The Autolev commands that fol- 
low are used to generate expressions for the nonlinear velocities of the centers of mass for the pushrod, 
arm, and flotor; i.e., IS PR, through the first actuator: 


V_° a _ Sa> = 0> , 

( 37 ) 

V2PTS (S<a, P a ,C> a ,A2 a ) , 

( 38 ) 

V2PTS (S a , P a , O a , P a o) , 

( 39 ) 

V P3.0 N>=V P3.0 S a > , 

( 40 ) 

V2PTS (S a ,A a ,A2<a, Fl a ) , 

( 41 ) 

V Fl a N>=V Fl a S a > , 

( 42 ) 

V2PTS (S a , A a , A2 a , Aao) , 

( 43 ) 

V A a o N>=V A a o S a > , 

( 44 ) 

V2PTS (S a , F a , Fl a , Fo) , 

( 45 ) 


V Fo N>=V Fo S a > . 


( 46 ) 
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4.7 Accelerations of the Centers of Mass for the Rigid Bodies 

Taking the time derivatives of the respective linearized velocity vectors yields expressions 
for the linearized accelerations of the centers of mass for each rigid body. This is accomplished with 
the following Autolev commands for the first actuator and in a similar manner for the additional seven: 


1* 

00 

0 

1 

s 

V 
II 

o 

V 

(49) 

1* 

O 

s 

V 
II 

o 

V 

(50) 

A2PTS (N, Pa,Oa,A2a) , 

(51) 

A2PTS (N, Pa,Oa, Pao) , 

(52) 

A2PTS (N, Aa, A2a, Fla) , 

(53) 

A2PTS (N, Aa, A2a, Aao) , 

(54) 


and 

A2PTS (N, Fa, Fla, Fo) . (55) 

4.8 Partial Velocities and Partial Angular Velocities 

Partial velocities and partial angular velocities are required with Kane’s method for determining 
the generalized active and generalized inertial forces. These calculations can be made explicitly by use 
of the Autolev PARTIALS command for a check of hand calculations. However, in the development of 
the Autolev model, the calculations can be performed automatically without the PARTIALS command 
by using the commands FR ( ) and FRSTAR ( ) . These commands calculate the generalized active and 
generalized inertial forces, respectively. 


4.9 Generalized Active Forces 

As mentioned in section 4.8, the generalized active forces and generalized inertia forces 
are calculated automatically by the FR ( ) and FRSTAR ( ) commands. Onorbit; i.e., neglecting the 
effects of gravity, the system is acted on by forces and moments due to the Lorentz coil, umbilical, 
and direct disturbance. Of these loads, only the forces and moments due to the actuator springs, 
and direct disturbances are considered. 

4.10 Contributions Due to the Pushrods 

The contributing loads on each pushrod are moments generated at the pertinent upper and lower 
stingers: 
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TORQUE (Pa, 0 . 001*Q2*Pal>+0 . 001*Q3*Ta3>) 


(56) 


and 


TORQUE (Pa, 0 . 001*Q4*Sa3>+0 . 001*Q5*Ra2>+0 . 001*Q6*Qal>) . (57) 

4.11 Contributions Due to the Actuator Arms 

The forces and moments acting on the first actuator arm are due to the respective Lorentz 
coil, located at A3 a; the flotor, through the first cross flexure; and the respective pushrod, through the 
upper stinger. Although the control input forces and moments are also contributing loads, they are not 
considered for the present Autolev models. 

The remaining contributing loads are due to the flotor and pushrods. The reaction forces between 
these elements are noncontributing. The contributing loads are 

TORQUE (Aa, 0 . 001*Ql*Fa3>) (58) 


and 


TORQUE (Aa, 0 . 001*Q2*Pal>+0 . 001*Q3*Ta3>) , (59) 

where the value, 0.001, denotes the cross-flexure stiffness (eq. (58)) and upper-stinger spring 
stiffness (eq. (59)) measured in lbf-in/rad. (Note the coupling between the control inputs and 
the generalized coordinates.) 


4.12 Contributions Due to the Flotor 

The actuators transmit loads to the flotor at the actuator cross flexures. These loads for the 
first actuators consist of a force applied at a point Fla and a moment exerted through the cross flex- 
ure. Since the force is noncontributing, it can be ignored in the analysis. The moment due to the cross- 
flexure spring is 


TORQUE (F, 0 . 001*Ql*Fa3>) , (60) 

where 0.001 is the cross-flexure spring stiffness measured in lbf-in/rad. The loads for the other actuators 
can be considered analogously. 

Unknown direct disturbance forces also act on the flotor. These can be represented as 
a force-moment pair, where the force is assumed to act through the flotor mass center, Fo. Let 
FORCE (Fo, ...) and TORQUE (F, ...) represent, together, this force-moment pair, where the ellipses 
indicate the magnitude and direction information, which must be supplied for an actual simulation 
run. Then the generalized-active-force contribution due to this pair could be determined explicitly, 
if desired, by the command 

GAF ( ) =DOT (V_Fo_N>, FORCE (Fo, . . . ) ) +DOT (W_F_N>, TORQUE (F, . . . ) ) . (61) 
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4.13 Generalized Inertial Forces 


As with the generalized active forces, the generalized inertial forces are also calculated auto- 
matically by Autolev in response to the FRSTAR ( ) command. Autolev uses the previously-calculated 
angular velocities and angular accelerations along with the user-defined inertial properties for these 
calculations. 

The central principal moments of inertia of the first actuator arm are defined as follows 
(using the designation Aa to indicate actuator number 1): 


INERTIA Aa,7e-5,7e-5,7e-5 . 


(62) 


Similarly, for the first pushrod: 


INERTIA Pa, 7e-4, 7e-4, 7e-4 . (63) 

The corresponding inertias for the remaining actuator arms and pushrods are defined analo- 
gously. These values for the actuator moments of inertia are reasonable approximations, and they can 
be modified as better data become available. The inertias for the ISPR and the ISS are entered respec- 
tively as 


INERTIA F, 66. 8, 166. 8, 71. 5, -0.03, 0.02, 1.17 , 


(64) 


and 


INERTIA S, 6875, 36875, 36875, -6, 4, 234 . (65) 

The centers of mass for the components are as entered in equations (33), (35), and (36). Again, 
once the mass properties have been defined, the generalized active forces are automatically calculated 
by Autolev in the course of evaluating FRSTAR ( ) . 

4.14 Constraint Equations 

All together, there are a total of 54 constraint equations: 48 kinematical and 6 dynamical. Since 
the complete set of equations for a system with n generalized coordinates consists of 2 n equations and 
since the ARIS system has 48 such coordinates, 42 more equations are needed to complete the set. These 
additional equations are the holonomic constraint equations (in nonholonomic form) for the dependent 
generalized speeds, Ui , ( i = 1 , ... , 4 8 ) . 

Since the velocity of the flotor center of mass, Fo, is the same irrespective of the actuator path 
chosen for describing its position, a set of constraint equations can be written in vector form for the first 
two actuators: 


DOT (V Fo Sa>, Fl>) -DOT (V Fo Sb,Fl>)=0 , 


( 66 ) 
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(67) 


DOT (V Fo Sa>, F2>) -DOT (V Fo Sb,F2>)=0 , 


and 


DOT (V_Fo_Sa>, F3>) -DOT (V_Fo_Sb, F3>) =0 . (68) 

There are an additional 3 equations for each of the remaining 6 actuators for a total of 21 con- 
straint equations; e.g., for the third actuator, one would substitute the lowercase “c” for “b” in 
equations (66)-(68). 

Since the angular velocity of the flotor, F, is the same, irrespective of the actuator path cho- 
sen, another set of constraint equations can be written in vector form for the first two actuators: 

DOT (W_F_Sa>, Fl>) -DOT (W_F_Sb>, Fl>) =0 , (69) 

DOT (W_F_Sa>, F2>) -DOT (W_F_Sb>, F2>) =0 , (70) 

and 

DOT (W_F_Sa>, F3>) -DOT (W_F_Sb>, F3>) =0 . (71) 

Similarly, one can write equations that compare the flotor angular velocity through the first 
actuator with the flotor angular velocity through the third actuator (substituting “c” for “b” as 
above), and so forth. This process yields the final 21 constraint equations. To finish implementing 
these 42 relationships, Autolev uses the command 

CONSTRAIN (DEPENDENT [U7 , U8 , U9 , U10 , Ull , U12 , U13 , U14 , U15 , U1 6 , 

U17,U18,U19,U20,U21,U22,U23,U24,U25, 

U2 6 , U2 7 , U2 8 , U2 9 , U3 0 , U3 1 , U32 , U33 , U34 , 

U35 , U3 6 , U37 , U3 8 , U3 9 , U4 0 , U4 1 , U42 , U4 3 , 
U44,U45,U46,U47,U48] ) , (72) 

After this final command has been executed, the 42 dependent generalized speeds are no longer 
used in the Autolev model. 
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5. COMPARISON BETWEEN NONLINEAR AND TRUTH MODELS 


To compare the truth model and the nonlinear model, the cross-flexure stiffness and stinger- 
spring stiffness were set to negligible levels, 0.001 lbf-ft/rad (0.00021 lbf-in/deg). Various test inputs 
were applied at selected target points to compare these models. Two such test inputs were (1) a constant, 
1-lb force with a constant direction relative to the IS PR, and (2) an oscillating force with an amplitude 
of 1 lb and an angular frequency of 0.2 rad/s. For the tests presented below, the target point was arbi- 
trarily displaced 6 in per direction, Fl>, F2>, and F3>, from the flotor center of mass. The direction 
of application for each force is along the Fl> direction; i.e., along the ISPR x-axis. 

The first comparison between the nonlinear and truth models (fig. 1) shows the response of each 
model to an oscillating input force applied directly to the ISPR. The respective model behaviors are very 
similar throughout the 10-s simulation. The differences between the two can be attributed to round-off 
errors in the fixed rotations between the ISS, ISPR and actuator-fixed coordinate systems, and some 
minimal effect due to the actuator springs. In this comparison, the actuators are not fully extended, 
so there are no noticeable effects from the constraint equations. 



-♦—Nonlinear: x-Axis Nonlinear: y-Axis ▲ Nonlinear: z-Axis 

Truth: x-Axis x Truth: y-Axis — Truth: z-Axis 

Figure 1. Comparison No. 1 between truth and nonlinear models. 


The second comparison (fig. 2) shows the response of each model to a constant input force 
applied to the ISPR. In this model, the response of the nonlinear model tracks the response of the truth 
model for =4.4 s, or until the mass center has moved «3 in (0.25 ft). At this point, an actuator linkage 
is fully extended, which produces a jerk, and the motion reverses direction. 
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Figure 2. Comparison No. 2 between truth and nonlinear models. 



This explanation for the behavior was further verified by the visualization afforded by the 
ENVISION model. At the time in the simulation when the ISPR center of gravity reversed direction; 
i.e., 4.4 s, the actuator configuration predicted by Autolev was recorded. That data were then entered 
into the ENVISION model, which clearly showed an actuator arm and pushrod having just reached full 
extension (fig. 3). 




Figure 3. ENVISION visualization of ARIS actuators near 
full extension. 
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6. LINEARIZED ACTIVE RACK ISOLATION SYSTEM MODEL 


With the full nonlinear model of ARIS completed, linearization is relatively easy using the 
Autolev system. The automatic linearization of the generalized active and generalized inertial forces 
is accomplished with the AUTOTAYLOR command. The arguments for AUTOTAYLOR describe 
the order of the terms of a Taylor-series expansion, here the zero and first order terms and the argu- 
ments about which the Taylor-series expansion is performed. For the first 10 generalized coordinates 
and generalized speeds, 

AUTOTAYLOR ( 0 : 1 , Q1 , Q2 , Q3 , Q4 , Q5 , Q6 , Q7 , Q8 , Q9 , Q1 0 ) , 

AUTOTAYLOR (0:1,01' , Q2 ' , Q3 ' , Q4 ' , Q5 ' , Q6 ' , Q7 ' , Q8 ' , Q9 ' , Q1 0 ' ) 
AUTOTAYLOR (0:1,U1,U2,U3,U4,U5,U6,U7,U8,U9,U10) , 
and 

AUTOTAYLOR (0:1, Ul' , U2 ' , U3 ' , U4 ' , U5 ' , U6 ' , U7 ' , U8 ' , U9 ' , U10 ' ) . (76) 

This is continued for the remaining 38 generalized speeds and generalized coordinates. The result 
is a linearized form of the ARIS dynamics. 


(73) 

(74) 

(75) 
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7. COMPARISON BETWEEN LINEARIZED AND NONLINEAR MODELS 


The small-angle theorem is based on a truncated Maclaurin- series expansion of the expressions 
for sine and cosine (about 0=0): 


and 


Q 3 

sin 6 = 0 

3 



0 2 0 4 

COS0 = 1 1 ... . 

2 4 


(77) 


(78) 


Depending on the application, this approximation is accurate for values of 6 less than =17°. 
In the present application, the constraint equations are most affected by the linearization, and values 
of 6 less than =8° produce small-angle approximations that remain valid. 

To ensure that the ISPR remains centered in its rattle space, and that the small-angle approxi- 
mations remain valid, the spring stiffness at the cross flexure and upper and lower stingers were next 
increased to nominal levels. The spring rates are set to 1 lbf-ft/rad (0.21 lbf-in/deg). In figure 4, one 
can see how the effect of the increased actuator spring constants changes the ISPR response to the 
applied force. The ISPR motion in the case for the moderate actuator spring constants is relatively 
small compared to the ISPR motion when the actuator spring constants are negligible. 



— Low: x-Axis ■ Low: y-Axis -a- Low: z-Axis 

Moderate: x-Axis Moderate: y-Axis • Moderate: z-Axis 

Figure 4. Low versus moderate actuator spring stiffness— 
response to constant test input. 
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Figure 5 shows the response measured in actuator number 1 angles when the spring constants are 
negligible and the ISPR is subjected to a constant input force. At =4.4 s, the actuator angles exceed the 
region bounded by the ±0.2 rad, where the constraint equations are valid. This point is further empha- 
sized in figure 6, which shows that the linearized model has very good correspondence with the nonlin- 
ear model until the simulation reaches 4.4 s. After that, the linearized model fails because the constraint 
equations no longer hold. 



Figure 5. Actuator No. 1 angles relative to small-angle 
approximation (±0.2 rad)— low actuator 
spring rates. 



— ♦ — Nonlinear: x-Axis — Nonlinear: y-Axis — ^ — Nonlinear: z-Axis 

Linear: x-Axis — ae- Linear: y-Axis — Linear: z-Axis 

Figure 6. Linearized versus nonlinear model, low stiffness— 
response to direct, constant test input. 
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Figure 7 shows the response measured in actuator number 1 angles when the spring constants are 
increased to moderate levels (0.21 lbf-in/deg) and the ISPR is again subjected to a constant input force. 
In this case, the angles never exceed the region defined by ±0.2 rad. In fact, the angles for actuator 
number 1 never exceed ±0.1 rad. Figure 8 shows a comparison between the ISPR motion in the linear- 
ized and nonlinear models when the actuator spring constants are increased to moderate levels. There is 
a nearly exact match between these models. 



Time (s) 


Figure 7. Actuator No. 1 angles relative to small-angle 
approximation (±0.2 rad)— moderate 
actuator spring rates. 
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Figure 8. Linearized versus nonlinear model, moderate 

stiffness— response to direct, constant test input. 
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This important result shows that, for small angles, the linearized model very closely tracks the 
nonlinear model, which, in turn, closely tracks the truth model. With the actual hardware, the ISPR must 
remain near the center of its rattle space due to the snubber system (fig. 3). Further, the control system 
for ARIS is designed to keep the ISPR from even contacting the snubbers. Assuming that truth model, 
which is relatively simple, is correct, then the linearized model presented in part one is accurate within 
operational or rattle- space limits. 
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8. LINEARIZED REDUCED-MASS MODEL OF THE ACTIVE RACK 

ISOLOATION SYSTEM 


A further simplification can be made by setting the mass properties (masses and inertias) of the 
actuator arms and pushrods to zero, while retaining those of the ISPR. The arms and pushrods represent 
16 of the 17 bodies in the generalized inertia forces of the nonlinear and linearized models, but only 
about 1.5 percent of the mass of an ARIS-outfitted ISPR. By setting the mass properties of the arms 
and pushrods to zero, the number of expressions in the set of holonomic generalized inertial forces 
eq. (A-180) of part one of this TM series) is reduced by 288 (r =1,..., 48, It was shown 

that dependent and independent generalized speeds, u r and u s , respectively, related as 

6 

u r = ^ \ s u s (for r = 7,... ,6) . (79) 

s=l 

The nonholonomic and holonomic generalized inertia forces, F* and F*, respectively, relate as 

48 

F* = F* + X F*A rs (r = 1,...,6) . (80) 

s= 7 

Since Kane’s method yields only as many dynamical equations as independent generalized 
speeds, removing these 288 expressions from the holonomic generalized inertial forces produces 
a net reduction of 1728 (6x288) terms in the final system equations. 

Despite the great reduction in the number of terms in the generalized inertia forces, the sim- 
plified or reduced-mass dynamical equations maintain the integrity of the generalized active forces, 
including control and disturbance inputs, and the constraint equations. The final dynamical equations, 
therefore, are the same as those reported previously in part one, excluding the contributions to the 
generalized inertia forces from the arms and pushrods. 1 
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9. COMPARISON BETWEEN LINEARIZED AND LINEARIZED REDUCED-MASS MODELS 


Since the generalized inertia forces are the most complex expressions in the dynamical equation 
development, the omission of the 288 expressions for the arms and pushrods represents a substantial 
reduction in the overall size of the dynamical equations. 6 This simplification is made by simply setting 
the appropriate masses, moments, and products of inertia to zero in the Autolev code, and then recompil- 
ing the generated C program. 

Figure 9 shows the response measured, in linear displacement of the ISPR center of mass, when 
the spring constants are increased to moderate levels (0.21 lbf-in/deg) and the ISPR is again subjected to 
a constant input force. From this plot, it is apparent that neglecting the mass properties of the arms and 
pushrods have little effect on the model’s description of system behavior. This is as expected, for two 
reasons: (1) As noted previously, the masses of the arms and pushrods are relatively small in comparison 
to the overall mass of an ARIS-outfitted ISPR, and (2) Because the ARIS arms and pushrods are more 
closely linked to the (relatively fixed) space station than to the ISPR, the amplitudes of their displace- 
ments, velocities, and accelerations will be correspondingly smaller than those of the ISPR. 



— t— Nonlinear: x-Axis — B — Nonlinear: y-Axis a Nonlinear: z-Axis 

Linearized Reduced Mass: x-Axis — • — Linearized Reduced Mass: z-Axis 

X Linearized Reduced Mass: y-Axis 

Figure 9. Linearized reduced mass versus nonlinear model, 
moderate stiffness— response to direct, constant 
test input. 
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10. CONCLUDING REMARKS 


Part one of this series presented a high-fidelity, linearized, analytical model of ARIS suitable 
for controller design . 1 In this TM, part two, two COTS software packages were used to develop compar- 
isons among nonlinear, linearized, and linearized reduced-mass models of ARIS (all without umbilicals). 
Direct and indirect disturbances were applied to these three models, and the responses were compared. 

The linearized reduced-mass model presented in this TM captures the significant system 
rigid-body dynamics of the high-fidelity, nonlinear model of ARIS. However, being simpler, it is more 
suitable for controller design while using standard controller-design tools. The relatively low level of 
model uncertainty could translate, in principle, into a high level of derived-controller performance. 

The techniques employed in this work could be used for solving and visualizing a variety of 
complex, over-constrained dynamics problems that model the interactions among several rigid bod- 
ies. However, these techniques are particularly advantageous for microgravity studies, where testing 
a microgravity isolation system in a standard Earth-gravity environment is difficult, if not impossible. 
The models developed here can be used as virtual test hardware for new configurations of umbilicals 
and payloads. 

As the ISS nears completion, and a full-time microgravity research program develops, it 
is reasonable to expect that there will be an increased emphasis on the quality of the ISS micrograv- 
ity environment. Suitable high-fidelity models of ARIS permit the development of correspondingly 
high-performance controllers. Improved controllers, in turn, facilitate the provision of an improved 
microgravity environment for the ISS science community. 

Part three of this series presents the incorporation of ARIS umbilicals into the ARIS model, 
which is the final step in the dynamical modeling process. The completed nonlinear model could 
be used for simulation studies or, in linearized form, for controller design using MATLAB® or other 
popular design tools. 
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